Dipolar interaction between two-dimensional magnetic particles 
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We determine the effective dipolar interaction between single domain two-dimensional ferromag- 
netic particles (islands or dots), taking into account their finite size. The first correction term 
decays as 1/ D , where D is the distance between particles. If the particles are arranged in a regular 
two-dimensional array and are magnetized in plane, we show that the correction term reinforces the 
antiferromagnetic character of the ground state in a square lattice, and the ferromagnetic one in a 
triangular lattice. We also determine the dipolar spin-wave spectrum and evaluate how the Curie 
temperature of an ensemble of magnetic particles scales with the parameters defining the particle 
array: height and size of each particle, and interparticle distance. Our results show that dipolar 
coupling between particles might induce ferromagnetic long range order at experimentally relevant 
temperatures. However, depending on the size of the particles, such a collective phenomenon may 
be disguised by superparamagnetism. 

PACS numbers: 75.75.+a, 75.30.Ds, 77.80.Bh, 75.20.-g 
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I. INTRODUCTION 

In this paper we are interested in two-dimensional mag- 
netic particles interacting through the long range dipolar 
forces. These particles have a two-dimensional character 
in two respects: first, because they are platelet shaped, 
that is to say their thickness t is much smaller than their 
linear size L; second, because they are arranged on a two 
dimensional substrate. 

They can be obtained, e.g., growing by Molecular 
Beam Epitaxy (MBE) a magnetic element on a high sym- 
metry substrate. In this case, growth is driven by surface 
diffusion, nucleation and aggregation:!!! in the submono- 
layer regime, the magnetic overlayer is made up of an en- 
semble of atomically thick islands that generally are not 
uniform in size neither arranged in a regular array. The 
distribution of islands may be regular if nucleation (and 
therefore island formation) takes place on.a reconstructed 
surface) 2 ] or on a network of dislocations!] Alternatively, 
particles can be produced via lithographic techniques :□ 
in this case, they are much bigger in size and their dis- 
tribution is generally uniform. 

In the following, small particles obtained by MBE 
growth in the submonolayer regime will also be called 
islands and large particles obtained by lithographic tech- 
niques will also be called dots. Particle is a generic term 
for both cases. 

Each particle is made up of a large number M of 
spins which interact ferromagnetically through the strong 
intra-particle exchange interaction. In an island, M w 
10 2 — 10 4 , while in a dot AT may be several orders of 
magnitude greater. 

Sufficiently small particles are expected to be in a sin- 
gle domain state, even if their actual magnetic state may 
depend on several factors: the shape of the particle, the 
strength of the anisotropies, the single crystal or poly- 
crystalline character of the particle, and so on. In this 



paper we are assuming that particles are in a single do- 
main state and have a crystalline structure. Within these 
hypotheses, the magnetic state of an isolated particle is 
fixed, first of all, by the balancing between dipolar in- 
teraction (which has an easy-plane effect) and possible 
anisotropies favouring the direction perpendicular to the 
plane (z direction, in the following). In the absence of 
quartic and higher order anisotropies a canted configura- 
tion is impossible and the resulting effect may be easy- 
axis or easy-plane only. 

The effect of dipolar interaction between spins belong- 
ing to the same (ultrathin) particle has been studied 
in a previous paperfl We showed that in-plane shape 
anisotropy is weak: this means that if the magnetization 
of the particle is within the film plane, its orientation is 
expected to be settled by the symmetry of the underlying 
lattice (through magnetocrystalline anisotropies) rather 
than by the shape of the particle (through the intra- 
particle dipolar coupling). Such a feature was indeed 
experimentally observed in MBE-grown Co on Cu(100) 
ultrathin particles.O 

In this paper we aim to study the interparticle dipolar 
interaction and to address the following questions: 

i) Assuming the particle to be in a single domain state, 
how is the dipolar interaction between particles modified 
by their finite size? 

ii) Assuming the particles to be arranged in a regular 
array, is the dipolar ferromagnetic (FM) state a stable 
configuration? 

Hi) If the ground state is ferromagnetic, is it stable 
at finite temperature and what is the value of the Curie 
temperature? 

The previous questions are addressed in Sections II, 
IV , respectively; in Section [V| the conclusions are 
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account the finite size of the particles: 



FIG. 1: Ii and I2 are two generic interacting particles. Each 
spin of the first particle, located in R\, interacts with each 
spin of the second particle, located in ife- Using reference 
frames centered in the centers of mass of the particles, the 
spatial positions of the two spins are fi , r*2 so that the distance 
between the two spins can be expressed as R12 = R2 — Ri = 
D + f; 2 — fi = D + f. 



II. EFFECTIVE DIPOLAR INTERACTION 
BETWEEN PARTICLES 



In the following we are going to consider two particles 
7i,/2 of any shape (see Fig. 1), with linear sizes L\,L2 
and thicknesses t\,t2- Each particle is a discrete collec- 
tion of spins and it is supposed to be in a single domain 
state: each spin of the two particles is indicated by si 
and s*2, respectively. The effective dipolar interaction be- 
tween the two particles is evaluated by taking the inter- 
action between a spin si located in R\ , a spin S2 located 
in i?2 and summing up on all them: 



Ri R2 



Si ■ S 2 (Si • i?l 2 )(s2 • -R12) 

7? 3 



R 5 



(i) 

where fl = g 2 p? B (g is the gyromagnetic factor and pb the 
Bohr magneton) and Co is the interplane distance in the z 
direction. In the previous expression, we have supposed 
that the thickness of each particle is much smaller than 
its linear size, t <C L. In this hypothesis, E dip is just 
linear in the numbers of atomic planes, ti/co and t2/co, 
and i?i2 are two dimensional vectors. 

As explained in Appendix [A|, the two quantities ap- 
pearing in square brackets in Eq. ([!]) can be expanded in 
the ratio r/D, where D is the distance between the cen- 
ters of mass of the particles and r = r*2 — fi (see Fig. 1). 
The effective dipolar interaction between two particles at 
distance D takes the approximate form: 



(2) 
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S 1 -S 2 o (5i-D)(^,.D) 
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9X12 Si ■ S2 
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D 5 
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(4'--p)(4'--p 





(3) 



(4) 



Each particle behaves as a single spin S — Afs, where 
J\f is the total number of spins in the particle. If L is 
the linear dimension of a particle and t its thickness, 
denoting by ao the in-plane atomic distance, one has 
N = 7(L/ao) 2 (t/co), where 7 is a geometric factor, de- 
pending on the shape of the particle and the lattice struc- 

(2) 

ture. The correction terms included in E dip decay with 
distance as 1/-D 5 , whilst the usual dipolar interaction 

decays as 1/D 3 . More precisely, E^l is a factor X12/ D 2 

smaller than E^} , where X12 = h(Zi +I2) is the semisum 
of the "moments of inertia" X of the two particles (see 
Appendix |a|) . 

In the continuum approximation, we have the following 
expressions. For a square particle of side L, X = L 2 /6; for 
a circular particle of radius p, X = p 2 /2; for a triangular 
(equilateral) particle of side L, X = L 2 /12. 

Finally, we would like to remark that the isotropic cou- 
pling term, proportional to (1S1 -£2), has different correc- 
tion terms according to the orientation of the spins: see 
Eq. (|J) . In other words, such a term is no more isotropic 
once the finite size of the particle is taken into account. 

In Fig. 2 we compare the exact dipolar coupling 
E dip (symbols) with the zero order approximation E d ® p 
(dashed lines) and with the second order approximation 
(E d ®p + E^p) (full lines), in two cases: i) the particles 
are magnetized perpendicularly to the plane, along the z 
axis (positive coupling energies), and ii) the particles are 
magnetized in plane along the x axis (negative coupling 
energies). We can see that the second order approxima- 
tion is fairly good except at very small distances: the 
smallest allowed distance between centers (without su- 
perposing the particles) is D min = V~2(L + a ) = 29.7a , 
where L = 20ao is the side of the square particle. The 
smallest value of D plotted in Fig. 2 is D = 30ao- 

In order to understand which configurations are en- 
ergetically favoured, let us start by considering just a 
couple of spins Si, S2, located in plane along the x axis. 
The dipolar coupling can be generally written as 



Ey 



P (Q) 
E12 



p(2) 
E\2 



VLq(S\ ■ S2 

-ri2(2Si S2 ' 



3S? flf) 

Si ■ S2 — 5Si S2) 



where E^ is the zero-order coupling and E d ^ p takes into 



where the explicit expressions of tto and ^2, depending 
on the distance between spins and on the size of the parti- 
cles, are irrelevant. We observe that Si, S2 may be either 
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FIG. 2: Dipolar energy between two square particles, one 
monolayer thick, of side L — 20ao on a square lattice (21 
spins per side). The centers of the two particles have the co- 
ordinates (0,0) and {D/y/2,D/y/2). The sides of the squares 
are parallel to the axis x, y. Positive and negative energies 
refer to spins parallel to z and x, respectively. The exact cal- 
culation E dip (symbols) is compared with E d ® p (dashed lines) 

and (E d ° ip + E d 2 ip ) ( mu lines). Inset: the case of spins parallel 
to x, in a log-log scale. 



veritable spins s*i , s 2 or they may represent the effective 
spins of two particles S-l, S2- in the former case one has 
O2 = 0, while in the latter (I2 is the correction due to 
the finite sizes of the particles. 

Both and Ejtf are the sum of competitive inter- 

(2) 

actions. For the sake of definitcncss, let us consider E\ 2 ■ 
The sum 25 f | : S'f + S\ ■ S2 favours an antiferromagnetic 
(AFM) alignment of the spins in the z direction, perpen- 
dicular to their joining vector. The term — 55f 5f favours 
a ferromagnetic (FM) alignment along the x axis. The 
energy of the former configuration is — 31^2 and the en- 
ergy of the latter one is — 4^2- As proved in Appendix^, 
where a more detailed discussion is given, the latter con- 
figuration is the ground state indeed. 

In conclusion, two spins interacting through the dipo- 
lar coupling minimize their energy by ordering ferromag- 
netically along the joining line. 

We are now going to discuss the more complex casg 
of a two dimensional lattice of spins. It is well knownQ 
that in the presence of a direct exchange interaction, the 
system is ferromagnetic and magnetized in the plane, be- 
cause of an easy- plane effect of E dip . If such an exchange 
interaction is absent, the easy-plane effect survives, but 
the actual configuration in the plane strongly depends on 
the lattice structure Jj It is useful to explain the origin of 
such a dependence on the spin arrangement. 

We have seen that spins would like to point along the 
line joining them: in a two dimensional lattice it is impos- 
sible, of course, to fulfill this requirement for all couples 
of spins. It is possible, however, for a chain of spins: so, 
we can start by addressing the nature of the coupling be- 



tween chainsS0 Let us consider a ferromagnetic chain of 
spins along the y axis of the plane and evaluate the dipo- 
lar field H dip generated at a point at distance d. Spins 
are oriented in the +y direction and _ff dip = H dip y. In the 
continuum approximation, 



H dip — 



- %nsi 



-Hoc. 



■1 1 ^ 



+00 



■n)- 

1 1 



3+n 



to 



where r — \/y 2 + d 2 . 

Both integrals have the form J dy{^p^ — (3 - 
It is sufficient to integrate by parts the term 1/r 
prove that the integral vanishes for any n and for any d. 
Therefore, in the continuum approximation H dip = 0. 

An exact calculation on a discrete lattice gives a finite 
value for H dip , but its sign depends on the actual lattice 
structure. Therefore, the reason for the sensitivity of the 
ground state on the spin arrangement is clear. Dipolar in- 
teraction favours the formation of spin chains magnetized 
ferromagnetically along the chain: these chains are very 
weakly coupled, in the continuum approximation being 
even uncoupled. The sign of the coupling and therefore 
the nature of the ground state do depend on the lattice 
structure. 

In the next Section we are going to analyze the nature 
of the ground state and to study the spin-wave spectrum 
with respect to a ferromagnetic alignment. 



III. GROUND STATES AND SPIN- WAVE 
SPECTRA IN TWO-DIMENSIONAL LATTICES 

In the following we consider an orthorombic lattice (see 
Fig. 3a) with primitive vectors 



ui = Ma,-b), u 2 = h{a,b) 



(5) 



where a = 2D sm(<p/2), b = 2D n cos(c/>/2), and D is the 
rhombus side. For <j> = ^ and (f> = ^, the orthorombic 
lattice reduces to a triangular and to a square lattice, 
respectively. The angle <f> can be supposed smaller than 
^, because if <f> > § there is just an interchange between 
a and b (see Fjg_3a). 

It is knownElE2l that the ground state is ferromagnetic 
for 4> < (f) c ss 80° and antiferromagnetic for larger values 
of the angle <j>. It is noteworthy that we find c/> c to be 
almost unaffected when E^ p is added to E^ p : as a first 
approximation, the finite size of the particle has therefore 

negligible effects on the ground state. We can say as well 

(2) 

that the finite size of the particles (i.e. the term E dip ) 
reinforces the FM character of the ground state for a 
triangular lattice, tj> = j, and the AFM one for a square 
lattice, (j) = This statement will be substantiated by 
the analysis of the spin-wave spectra at the end of this 
Section. 
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(a), 



a 



of mass are located in Di and Dj, respectively. In the 
following we are assuming to have an array of identical 
particles, so that Xj = X. 

Taking the y direction as quantization axis, we per- 
form the Holstein-Primakoff transformation from spin to 
boson operators 



s, t 



S 

a 3 ), S* =S- a]a v 5? = \j -(a 3 + a)) 

Next, exploiting the translational invariance in the film 
plane, we introduce the Fourier transform 





ik-Dij 




k-Dij 



(8) 

where k = (k x ,k y ) is the two-dimensional in-plane 
wavevector ranging over the first Brillouin zone, gener- 
ated by the primitive vectors (see Fig. 3b) 



«i = M~,-|). « 2 = Mii) 



(9) 



and Nu is the total number of spins in the two- 
dimensional lattice. The spin-wave Hamiltonian takes 
the form 



FIG. 3: (a) Orthorombic lattice with coordinate axes directed 
along the rhombus diagonals, a and b; <f> is the rhombic angle, 
(b) Reciprocal lattice of the orthorombic lattice. The angle 
8k denotes the orientation of the in-plane wavevector k with 
respect to the x axis. 



Let us now calculate the frequency of spin-wave ex- 
citations with respect to a ferromagnetic configuration 
with magnetization directed in plane. We remind that, 
depending on the value of the rhombic angle (j>, the 
ferromagnetic ground state is different^ for (j> < 7r/3, 
the dipoles are oriented along the short diagonal of the 
rhombus (i.e. along the x axis in Fig. 3a) while for 
tt/3 < 4> < 4> c , along the long diagonal (y axis). 

The effective dipolar interaction between particles of 
spins whose centers of mass are located on the sites D 

of an orthorombic lattice is H. dip — Ti. d °n 
obtained from Eqs. (3,4) summing over all sites 



it is 



^EE 



Di Di 



Si ■ Sj ^(5j ■ Dij)(Sj ■ Dij) 



D 3 



D- 



H 



(2) 



Di Dj 



5/ • S 



D 5 



3X ? > 



Si ■ Sj (Si ■ Dij ) (Sj ■ ) 



D 5 



(6) 



where D^ = Dj — Di andX^ = is the semisum 

of the "moments of inertia" of two particles whose centers 



p = E A k a \H + \ B kiH a -k + 4 a %) 



The coefficients A? = A- 

k k 



(0) 



A { ? ] and B % 



(10) 



can be expressed (see Appendix C) through the dipolar 
sums (a, j3 = x,y, z) 



v {n hk) = y -4— 



D 

Di y 



l-(3 + n)- 



DZD 



2i v_ 

D 2 



Jk-Di 



(11) 



where for n = and n = 2 one has respectively the zero 
and second order expressions in the ratio L/Dq between 
the linear dimension of the particle and the interparticle 
distance. The spin-wave energy is 



{ A k- B k){ A k + B k 



1/2 



(12) 



In the continuum limit k — > 0, one finds the approximate 
analytic expression 



ns 



(A 1 -- 

V a 



b 

/ Auk \ 1 1/2 

(a (,/>) + — cos^fe + O^ 2 ))] (13) 



0(k 2 



where Ao and Ai, defined in Appendix C, depend on 
the angle 4>. In particular, A vanishes for the highly 
symmetric cases <p — tt/3 (triangular lattice) and 4> = n/2 
(square lattice). The angle Ok defines the orientation of 
the in-plane wavevector k with respect to the x axis. 
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It is worth stressing that eg has the same structure 

both at zero and second order, because H^J does not 
contribute to the linear term in k. In particular, we have 
that the linear term in the quantity (A^ + -Bg) vanishes 

when k is oriented along the magnetization (Ok = 7r/2). 
Its sign is therefore decided by the quadratic term, which 
is different for different lattices. It is positive for the tri- 
angular lattice (see Fig. 4) and it may be either positive 
(Fig. 5b) or negative (Fig. 5a) for the square lattice, de- 
pending on the orientation of the magnetization. 

In Figs. 4,5 we show the dispersion relation of the spin 
waves, as obtained from Eq. (12), in the case of a trian- 
gular lattice (</> = 7r/3) and of a square one (4> — 7r/2), 
for different orientations 9 k- 

In the remaining part of this Section we are making 
general comments on these results. 

Dipolar interaction, since it couples the spins to the un- 
derlying lattice, is not rotationally invariant, contrary to, 
e.g., the Heisenberg interaction (TL = —Jsi ■ sa). There- 
fore, the energy of a generic spin configuration depends 
on how the spins are oriented with respect to the lattice. 
However, lattices with four-fold (4> — 7r/2) or six- fold 
(<f> = 7r/3) symmetries are special, in the sense that the 
energy of certain configurations (the ferromagnetic one, 
for example) is rotationally invariant. This means that 
the dipolar ferromagnetic ground state of a triangular 
lattice has a continuous degeneracy. 

The coupling between spins and crystal lattice man- 
ifests itself in two ways (we are considering the highly 
symmetric cases, Ao(0) = 0). First, the spin-wave en- 
ergy eg depends on the orientation of the wavevector in 

the continuum limit k — > as well. According to Eq. (13), 
eg ~ | cos 6*fc Ifc 1 / 2 , showing that the (^-dependence of eg 

is maintained in the limit k — > 0. 

Second, even if the energy of the FM state is degenerate 
with respect to its orientation in the plane, the spin- wave 
energy is not (see Fig. 5). This is true for the triangu- 
lar lattice as well (not shown), but it has most striking 
consequences for the square lattice. In this case, the FM 
configuration is known to have a higher energy than the 
structure where spins are ferromagnetically coupled along 
lines and antiferromagnetically coupled between neigh- 
bouring lines. □ The FM state, if oriented along an axis of 
the square lattice (the typical configuration considered in 
the literature) is not even locally stable. This is clearly 
shown in Fig. 5a where the square of the spin-wave en- 
ergy is plotted, in the case of k and M parallel to the 
a side of the square lattice. The square of the energy 
is negative, signaling that the corresponding FM state is 
unstable. However, if the magnetization is oriented along 
the diagonal, the FM state is no more unstable, as shown 
in Fig. 5b (we show just the case of k parallel to M, but 
no instability appears for any value of 9 k). 

(2) 

Figures 4 and 5 display the effect of on the spin- 
wave spectrum: the energy increases for any value and 
orientation of k. This result corroborates our statement 
that the finite size of the particles reinforces ferromag- 
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FIG. 4: Spin-wave dispersion curves calculated for the tri- 
angular lattice (rhombic angle <f> — 7r/3) both numerically 
(Eq. symbols) and in the continuum limit (Eq. (|l3|), 

lines). Solid circles and full lines refer to H^; open circles 

and dashed lines refer to (H^ + H^), with X = 0.1 (units 
with Do — 1 are used). The magnetization is assumed to lie 
along the y axis and three different propagation directions are 
reported: (a) 6 k = 0; (b) 8 k = tt/6; (c) 6 k = tt/2. Note the 
different periodicities of the spin-wave energy. 



netism in the triangular lattice (Fig. 4). For the square 
lattice we have showed that ferromagnetism is locally sta- 
ble if M is oriented along a diagonal of the square lattice 
(Fig. 5b) and locally unstable if M is parallel to a side of 

the square lattice (Fig. 5a): E^J reinforces the stability 
in the former case and the instability in the latter one. 
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k y b 



W elc = nS 2 /Dl « {gji B s) 2 {L/a )y/cof/Dl the effec- 
tive dipolar interaction between different particles. 



The convergence of the integral (14) in k = is easily 
proved .□ For generic orthorombic lattices, it is guaranteed 
by the gap in the dispersion curve, eo 7^ 0. For square 
and triangular lattices, the spin- wave energy vanishes for 
k — > and one can expand the exponential on the de- 
nominator of Eq. ( |l4| ) , because the possible divergence is 
infrared-like. Thus, the integral giving the temperature 
dependent spin deviation is found to converge, provided 
that the positivgt3 0(k 2 ) terms in the spin- wave energy 
are taken into account. Eq. ([u]) can be rewritten as fol- 
lows: 




k y b 



FIG. 5: The same as in Fig. 4, but for the square lat- 
tice (rhombic angle <f) = 7r/2). (a) When magnetization and 
wavevector k are directed along a side of the square lattice, 
the square of the spin-wave energy is found to be negative, 
signaling the instability of the ferromagnetic configuration, 
(b) When magnetization and wavevector k are directed along 
the diagonal of the square, the ferromagnetic state is found 
to be metastable. 



IV. THE CURIE TEMPERATURE 



ss 
"5 



T 

W eU jit 



dO 



qdq 



ciq cos 2 9 + C2q 2 



(15) 



where q = kDo, 9 — 9k, c\ and ci are positive constants, 
and qu ~ 1- The double integral can be easily evaluated, 
giving 



Tr 



W„ s ■ — 

Z7T 



hi 



<7m + \Jq\i + (ci/c 2 ) 2 



-I -1 



ci/c 2 



(16) 



In the limit ci = 0, 



C2<7 2 reproduces the disper- 



sion curve of the Heisenberg ferromagnet and the Curie 
temperature Tc "vanishes, in agreement with the Mermin- 
Wagner theorem. 113 

Finally, it is interesting to discuss the case where a 
uniaxial single-ion anisotropy A, favouring the in-plane y 
axis, is present in the system. The effect of the finite size 
of the particle is straightforward, in this case: 



y\2 



(17) 



D, 



The Curie transition temperature of an orthorombic 
lattice of particles with magnetization directed in plane 
along the y axis can be estimated in the framework of 
spin-wave theory. The relative deviation of the magneti- 
zation from the saturation value takes the formQ 



5S 
5" 



N n S ^ K k k/ 



V 2 



{2ir) 2 S 



2eg 



(14) 



where V2 is the volume of the two-dimensional unit cell 
and the integration is over the first Brillouin zone. The 
first term on the r.h.s. of Eq. (14) gives the temperature 
dependence of the magnetization, while the second term 
represents the zero-point spin deviation, which can be 
safely neglected. A rough estimate of the Curie temper- 
ature T c is obtained by imposing that SS(T C )/S « 1. 

Since the only energy scale in the problem is given by 
the dipolar interaction, Tc is expected to be of order 



where A = \/Af and, as usual, S = Afs. 

The effect of E ani on the spin- wave energy is simply to 
add the quantity 2A to A^; for large particles, this con- 
stant factor dominates the dipolar terms in Eq. (13). Be- 
cause of that, spin- wave approximation is no more suit- 
able and Eq. (|14| ) can not be used to evaluate Tc'. in fact, 
it would give a Curie temperature Tc ~ 2J\f\, which di- 
verges when A — > 00. 

On physical grounds, we expect that Tc is an increas- 
ing function of A, but in the limit of strong anisotropy, Tc 
is always of order W eB . An analogy can be done with the 
three dimensional Heisenberg model in the presence of an 
easy-axis anisotropy A: Tc increases with A, but Tc(0) 
and Tc (00) are of the same order of magnitude and are 
both of order J, the exchange coupling constant. In our 
case, W„ft replaces J and the dimension of the system 
is two instead of three (the two-dimensional Heisenberg 
model is not ordered at finite temperature in the absence 
of anisotropy) . 

In order to corroborate our argument, Tc will be cal- 
culated in the mean field approximation, starting from 
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the Hamiltonian 



H 



MF 



E 



Ms?) 



(18) 



by 



where £ = ^n(S«)V^(0). The magnetization is given 

(s v ) = ^ M =- s 



E S m=s Me-« M "^ M2 )/ T 



M=-S 



-({M-AAf 2 )/T 



(19) 



where M, M 2 denote the eigenstates of S y , (S y ) 2 respec- 
tively. As T -> T c one has -> so that £ -> and 
the exponential can be expanded. Thus (S y ) ra —-^R 
where i? is defined as 



i? 



M=-S 



e T c 



((S v ) 2 )t c - (20) 



For A = one has R = §«S(«S + 1), so that T C (A = 0) = 
fa[-V$>(p)]S(S + l)KW«. 

For A 7^ the mean field critical temperature is 



T C (X) = 3 
T c (A = 0) S(S + 1 



■((s y ) 2 ) 



Tc 



(21) 



where the mean on the r.h.s. must be calulated numeri- 
cally. Now we observe that in the limit A — > oo, one has 
((S y ) 2 ) t c -> S 2 , so that the ratio T c (oo)/T c (0) tends 
to the finite value 35/(5 + 1). We conclude that for any 
value of A, the Curie transition temperature remains of 
the order of the effective dipolar interaction W eff , with a 
prefactor changing by a factor three as A increases from 
to oo. 



V. DISCUSSION 

Let us get back to the three questions formulated at 
the end of the Introduction. 

i) What is the effective dipolar interaction between sin- 
gle domain particles? 

In the hypothesis that spins are strongly coupled fer- 
romagnetically inside each particle, it is straightforward 
to define an effective dipolar coupling between (micro- 
scopic) spins s of two particles at distance D: W et{ = 
Q(sAf) 2 /D 3 , where ft = (gfJs) 2 and M is the number of 
spins in each particle. If L and t are respectively its linear 
size and thickness, W off ra (gnBs) 2 {L / 'ao) 4 (t/ 'c ) 2 / 'D 3 . If 
a sipgle-ion anisotropy A is present, its effective value is 
justlil Acff = AAA ra \(L/a ) 2 {t/c ). 

The full dipolar interaction between particles can be 
expanded in (even) powers of (L/D) m . The first correc- 
tion term (m = 2) gives an interaction decaying as 1/D 5 
with the distance between particles. A couple of remarks 
are in order here. First, the 'purely' spin term in the 



dipolar interaction (i.e., the term not coupled to the lat- 
tice) is no more rotationally invariant: the z-components 
are more strongly coupled than in-plane components. 

n (2) 

Second, preserves two important features: it is mini- 
mized when spins are aligned ferromagnetically along the 
joining line, and the dipolar field generated by a 'contin- 
uum' line of spins aligned along the line in a point outside 
the line, vanishes. 

ii) What is the dipolar ground state of an ordered array 
of magnetic particles ? 

We have considered the class of orthorombic lattices, 
which comprises the triangular and the square lattices. 
It is known that in the case of a lattice of spins, the 
ground state for the six-fold symmetry is ferromagnetic 
and for the four-fold symmetry has zero net magnetiza- 
tion. These results are not modified when single-domain 

(2) 

particles replace single spins and is considered in 

addition to E^. The effect of is to reinforce, in 

some sense, the effect of E^. In particular, therefore, 
a triangular lattice of two dimensional particles interact- 
ing through the dipolar interaction has a ferromagnetic 
ground state. 

Hi) What is the finite temperature behaviour of a in- 
plane ferromagnetically ordered array of particles ? 

In two dimensional systems, long range order at finite 
temperature is not certain: however, the long range dipo- 

(2) 

lar forces are known to grant it. In this respect, E\^ is of 
minor importance, because it decays as 1/D 5 and conse- 
quently the dipolar sums in the fc-space do not contribute 
to the terms linear in k, but to the quadratic terms only. 

In the absence of anisotropies, it is elementary that the 
Curie temperature Tc is of order of the effective dipolar 
coupling, Tc » Weft, because it is the only energy scale 
in the problem. If in-plane easy-axis anisotropies are 
present, we have shown that Tc is expected to increase, 
but not to change in order of magnitude: according to 
mean field theory, Tc increases by a factor three pass- 
ing from the 'weak' anisotropy regime into the 'strong' 
anisotropy regime. 

Denoting by w = (gfiBs) 2 / 'ag the dipolar coupling 
between microscopic spins on a two-dimensional lat- 
tice with atomic distance do, the effective dipolar in- 
teraction between particles of linear dimension L and 
thickness t at distance Do can be rewritten as W eff — 
w(L/ao) 4 (t/co) 2 /(D /ao) 3 . Thus, the Curie tempera- 
ture Tc ~ Weff of an ensemble of magnetic particles 
may be significantly larger than the Curie temperature 
of a two-dimensional lattice of microscopic spins, which 
is of order w. However, L cannot be made larger than 
Do, because Do scales with L: at the best, therefore, 
T c » w(L/a )(t/co) 2 . 

Even if we have no definite evaluation of the numerical 
prefactor appearing in the previous estimate for Tc, we 
suggest that regular arrays of two-dimensional particles 
magnetized in plane might sustain long range order at 
experimentally relevant temperatures (see also Ref. fl4| ). 
However, this collective phenomenon may be masked by 
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the superparamagnetic behaviour of the single particle, 
appearing below the blocking temperature Tg rs Ae« ~ 
X(L/a ) 2 (t/c ). 

The condition Tg > 7c is equivalent to X/w > 
rjL 2 t/DQ, where 77 is an unknown numerical factor. If 
such a condition is satisfied, the Curie phase transition at 
T = Tc is not visible because thermodinamic equilibrium 
cannot be attained below Tg. So, if we decrease temper- 
ature from the high-T region each particle becomes su- 
perparamagnetic when dipolar forces are still unable to 
induce a long range order in the system. When T = Tc 
dipolar interaction comes into play, but the magnetiza- 
tion of each particle is frozen and the system is unable 
to attain equilibrium. In the opposite case, Tc > Tb, a 
phase transition at Tc should be visible. 

Recently, the system Co/Cu(001) has_drawn the at- 
tention because it has been suggestecOiij that for t < 
1.8ML (ML=monolayer) this system displays a dipolar 
induced ferromagnetic order. Tc(t) is seenEJ to be finite 
and increase from T C (1ML)=25K to T C (1.8ML)=200K, 
where it has a sudden jump, attributed to the percolation 
in the second layer. 



The growth morphology o: 
complicated by the alloyin 



system Co/Cu(001) is 
between the two ele- 



ments, which mainly takes place in the first layer (75% 
in Co, 25% in Cu). However, with such a high percent- 
age of Cobalt, it is hard to suppose that the first layer 
is made up of an ensemble of disconnected Co-islands 
which interact only through long range forces (in Ref. [l6| 
the percolation threshold is theoretically estimated to be 
of order 60%). More likely, an infinite cluster of Cobalt 
does exist in the first layer and therefore the magnetic be- 



haviour of the system-. (and the value of Tc) follows from 
the combining effectt3 of the direct exchange interaction 
between spins and long range forces. 

Finally, we would like to mention an additional diffi- 
culty in the interpretation of experimental data concern- 
ing an array of magnetic islands: the random character 
of deposition and diffusion gives rise, to a non-uniform 
distribution of sizes and positions.EJ This fact, along 
with frustration due to dipolar interaction, makes diffi- 
cult even the determination of the ground state, because 
the system has a glassy behaviour with a lot of metastable 
states. 



APPENDIX A: MULTIPOLAR EXPANSION OF 
INTERPARTICLE INTERACTION 

We start from Eq. (]]]), assuming single monolayers .0 
Summations on i?i,2 are replaced by sums on ri,2 (see 
Fig. 1 for notations): 



E. 



.-to: 



S\ ■ S 2 (Si • #12) (s 2 ' R12) 

R 3 

rt 12 



R 5 

n 12 



We make the expansions: 



R 3 



1 

D 3 



f- D 3 r 2 15(r-L>) 5 
' D 2 ~ 2D 2 + ~2 D 4 



(Al) 



(A2) 



J 



(si ■ #12X5*2 R 



V2 



R 5 u 



1 



r ■ D 



+ [(Si • D){g 2 ■ f) + (s 2 ■ D)(Si -r)] 

I 



5 r 2 35 (r ■ D) 
2D 2 + 



_D 4 



+ (si ■ r)(s 2 ■ f) 



f ■ D 



(A3) 



The following expressions are easily calculated: 

1 = M E 1 = ^2 



EE^=^({rf) 



2(ri) ■ (r 2 » 



EE^ 1 ' ^2 ■ r) = MxM2[{(V x ■ ri)(F 2 • r x )) 
+ ((Vi ■ r 2 )(V2 ■ r 2 )) - (Vi ■ <r 2 ))(^2 • 

-{vi-{?i)){V2-{m 



In the previous expressions, A/1,2 are the number of 
spins in the particles I\,2 and V\ t 2 are generic vectors. 

It is always possible to choose the origin of the refer- 
ence system for a given particle in its center of mass, so 
that (fi) = and the "moment of inertia" li of a particle 
is Ti = (rf). Consequently, we have the following results 



2 U Xj 



{(D-n) 2 ) 
((Vi ■ n)(V2 ■ n)) = 



(A4) 
(A5) 



where V" is the in-plane component of the generic vector 
V. 

If we define 2i 2 = +2" 2 ) and Si = A/istj, we obtain 
the following expression for the effective dipolar interac- 
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tion between two particles at distance D: 



E dlp = 



2 D 3 
2 D 3 



91^ 
2D 2 

2D 2 



3 (S^-m-D) 
2 D 5 



5Xi2 
2 D 2 



(A6) 
(A7) 
(A8) 



APPENDIX B: MINIMIZATION OF THE 
DIPOLAR COUPLING BETWEEN TWO SPINS 

Let us consider two unitary spins S\ , S 2 located at a 
distance R± 2 along the in-plane a; axis, taken as the po- 
lar axis, while the z axis is perpendicular to the plane. 
Their orientations are defined by the polar and azimuthal 
angles 6i, ipi. 

The dipolar coupling can be generally written as 



E12 



p(o) 



E 



(2) 



~ rio(>s'i • s 2 — 3Sfs 2 ) 

+tt 2 {2S{S z 2 +S 1 -S 2 - 5SfS£) 

where f2 = \SIS 2 / R\ 2 and f2 2 = jlii^S 2 / ' R\ 2 (see 
Eqs. ||). 

Both terms E^ 2 and are minimized by the same 
configuration. First, let us treat the zero-order term. We 
have to minimize the function 

e[^ /Qq — sin 61 sin 9 2 cos{<pi—Lp 2 ) — 2cos6i cos 9 2 (Bl) 

By taking the derivatives with respect to ipi 2 we find 
that (<pi — 1P2) — 0, 7r or that one 6i at least must vanish. 
If, e.g., 61 = 0, it is straightforward to derive that 9 2 = 
as well. If both (9 12 are not vanishing, taking the deriva- 
tives with respect to them implies cos#i = cos #2 = 0, 
i.e. 6*1.2 = §• In simple words, it is sufficient to con- 
sider two kinds of configurations: i) the configuration 
where both spins are perpendicular to the joining vector 
and they are parallel or antiparallel; ii) the configuration 
where both spins are aligned along the joining vector. 
The former configuration, with antiparallel spins, corre- 
sponds to the minimization of (Si ■ S 2 ) and its energy is 
— fio- The latter configuration, with parallel spins, corre- 
sponds to the minimization of (—3SfS 2 ) and its energy is 
f2o~3fio = — 2r2 . We conclude that e[^ is minimized by 
the ferromagnetic configuration with spins aligned along 
their joining vector. 

(2) 

The minimization of E l2 proceeds along the same 
lines, with the minor difference that configurations with 
both spins perpendicular to the joining line are no more 
degenerate with respect to a global rotation around the 
x axis: the lowest energy one corresponds to antiparal- 
lel spins along the I direction and its energy is — 3r2 2 - 
On the other hand, the ferromagnetic configuration with 



both spins parallel to the x axis has the energy —ACl 2 , so 
the conclusion is unchanged. 



APPENDIX C: DIPOLAR SUMS FOR THE 
ORTHOROMBIC LATTICE 

For spins ferromagnetically oriented along the y axis, 
the coefficients of the spin- wave Hamiltonian in Eq. jlO| ) 
take the form 



4 



(0) 



ns 

k 

ns 



(CI) 



and 



AC 

k 



(2) 



k 

nsp: 



lvf}{k) + \v^(k)-V^{Q) + E^{k) 
\vf} (k) - (k) + iVf} (k) + £ W (fc) 



(C2) 

where Af is the number of spins in each particle and the 
dipolar sums are defined as (a, f3 = x, y, z) 



1 



Dj V 



3+n 



1- (3 + n)- 



D?,D 







D 2 



ik-Di 



£ {n \k) = E 



Ak-Di 



r\3+n 



(C3) 



We observe that in the ultr aflat particle limit, one has 
£(»)(&) » Di^ik) and D&Qe) » and the V aa 's can 
be approximately expressed 



Vi n j(k) » yW(£)-(2 + n)xW(2) 
V${k) « X(n)(k)-(2 + n)Y^(k) 

vif(k) m i w (fe) + yW(fe) 

in terms of the dipolar sums 



(C4) 



D 3 V 

(D y ) 2 



ik-Di 



ik-Di 



(C5) 



The latter sums can be numerically calculated in a very 
efficient way following ajnethod, developed some years 
ago by Benson and Mills,E3 similar to Ewald's one for the 
evaluation of lattice sums. For n — one has 

X^{k x ,k y ;a,b) = U^{k x ,k y ) + V^(k x ,k y ) (C6) 
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where U is a summation over the sites of the orthorombic 
lattice which are located at integer multiples of a and b 



and 



16 1 



3 b 



E ( 



m— — oo 



7rm 



V -) 2 K 7 . 



„ , . 7T771 k v 



while V refers to semi-integer multiples 



E 



7rm 
IT 



16 1 

y & 

fry 
2 



E« 

i=0 



„ l.,7rm 
2^+2)1- 



For n = 2 one has 

X^(k x ,k v ;a,b) = U^(k x ,k y ) + V^(k x ,k y ) (C7) 
where 

(2) 32 1 ^cos(^oi) 



/=i 



i 



E c 2 ? + ^) 3 ^ 



2a/|- + f| 



and 



' (fr'x i fry) 



, . . . +oo 

32 1 ^ cos 



15 ab 



E mhT + t) 3 ^ 



[k x a(l + |)] 

„ l.,7rm fc„ . 
2«G+2)|- + fl 



Here above, K2{z) and -^3(2) are the modified Bessel 
functions of second and third order, respectively. By 
symmetry reasons, one has 



A (2) (fc) » A (2) (0) + <3(fc 2 ) 

Y^ 2 \k) » y( 2 )(o) + o(fc 2 ) (cio) 

where for the triangular lattice one has Z(°'(0) = 
5.5170879/Dg, X<-^(0) = 3.3809493/Dj] and for the 
square lattice X^(0) = 4.5168109/Dg, X^{0) = 
2.5451291/£>q. From the previous expressions one 
obtains app roximate expansions for the quantities in 
Eqs. (Hill: 





a A + 


47rfc 


ns 


a 6 




a Ai- 


47rfc 


ns 


a 6 






27rfc 


ns 


a A 2 - 


ao 



^-^ sin 2 6» fe + 0(k 2 ) (Cll) 



where 

Ao = 
Ai = 
A 2 - 



r(°'(o)-i(°'(o) 



15, 



y( 2 )(o) ~x^(o) 



3y(°)(o) + |i 



7F (2) (0) + 2A (2) (0) 



2y(°)(o) - x(°)(o) 



4r^(o) -x( 2 )(o) 

(C12) 



Hence we observe that for a generic orthorombic lattice 
one has X^ffl) ^= FW(0), so that the dispersion curve 

1/2 ~* 

e fe = [( A k ~ B fc)(^fc + B k)] presents a gap for k -> 0. 
For the special cases of the triangular and the square 
lattice one has lw(0) = FW(0) by symmetry reasons; 
this implies that Aq = 0, so that a Goldstone mode is 



(k x ,k y ;a, b) = X ( - n) (k y , k x ; b, a) (C8) present in the dispersion curve: -> for fc -> 0. 



Finally, in the continuum limit k — > 0, one can obtain 
analytical expressions for the dipolar sums. Denoting by 
k the modulus of the two-dimensional wavevector k = 
(k x , k v ) and by 9k the angle that k forms with the x axis, 
we obtain 



ao 



yW(fe) » r w (o) 



(0), 



2tt 
ao 



l + icos(2^) 



1 - -cos(26> fe ) 
3 



(C9) 
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